…hen calling .to_array().sel(variable='band_1',

   drop=True), the resulting dask array has a different chunk structure than arrays created directly. The crosstab
  function failed because it pairs dask blocks positionally via to_delayed().ravel(), which requires aligned chunks.  To fix this:

  Added chunk validation for 2D inputs in crosstab, using the existing validate_arrays function (which is already
  used by zonal_stats). This function:
  1. Validates that shapes match
  2. Automatically rechunks values to match zones chunk structure when they differ